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Abstract 

We compare predictions of the mean-field theory of superconductivity 
for nearly-antiferromagnetic and nearly-ferromagnetic metals for cubic and 
tetragonal lattices. The calculations are based on the parameterisation of an 
effective interaction arising from the exchange of the magnetic fluctuations 
and assume that a single band is relevant for superconductivity. The results 
show that for comparable model parameters, the robustness of magnetic pair- 
ing increases gradually as one goes from a cubic structure to a more and more 
anisotropic tetragonal structure either on the border of antiferromagnetism 

or ferromagnetism. 
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I. INTRODUCTION 



One can expect that the effective interaction between quasiparticles in strongly corre- 
lated electron systems to be very complex. The interaction will depend obviously on the 
charge, but also more generally on the spin and current carried by the quasiparticles. On 
the border of long-range magnetic order it is plausible that the dominant interaction chan- 
nel is of magnetic origin and depends on the relative spin orientations of the interacting 
quasiparticles. 

It has been shown that this magnetic interaction treated at the mean field level can 
produce anomalous normal state properties and superconducting instabilities to anisotropic 
pairing states. It correctly predicted the symmetry of the Cooper state in the copper oxide 
superconductors and is consistent with spin-triplet p-wave pairing in superfluid ^He 
[for a recent review see, e.g., ref. 0]. One also gets the correct order of magnitude of 
the superconducting and superfluid transition temperature when the model parameters 
are inferred from experiments in the normal state of the above systems. There is growing 
evidence that the magnetic interaction model may be relevant to other materials on the 
border of magnetism. 

Thus far the magnetic interaction model has been explored in very simple cases. The 
most extensively investigated example is that of a nearly half-filled single-band in a square 
or cubic lattice. These studies have revealed a number of interesting features that are 
quite in contrast to those expected for conventional phonon mediated pairing. In the latter 
case, the interaction is local in space, but non-local in time, whereas on the border of 
magnetism, one expects the interaction to be strongly non-local in both space and time. 
For nearly antiferromagnetic metals the magnetic interaction is oscillatory in space and 
superconductivity depends on the ability of the electrons in a Cooper pair state to sample 
mainly the attractive regions of these oscillations. Because of the strong retardation in time, 
the relative wavefunction of the Cooper pair must be constructed from Bloch states with 
wavevectors close to the Fermi surface. Furthermore, the allowed symmetries of the Cooper 
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pair wavefunction are restricted by the crystal structure. The possibihty of constructing 
a Cooper pair state with maximum probabihty in the attractive regions of the magnetic 
interaction can be severely constrained by these requirements. Therefore, one expects that 
the robustness of magnetic pairing to be very sensitive to details of the electronic and lattice 
structures. 

On the border of ferromagnetism, one is not hampered by the oscillatory nature of the 
magnetic interaction which in the simplest model is attractive everywhere in space and time 
in the spin-triplet channel. At first sight this would seem to be the most favourable case for 
magnetically mediated superconductivity. However, the results of the numerical calculations 
presented in ref. indicate that the highest mean field for the cases considered is obtained 
for d-wave pairing in the nearly antiferromagnetic state in a quasi-2D tetragonal lattice. In 
this particular explained in ref. 0, it turns out to be possible to ideally match the 

Cooper pair state to the attractive regions of the magnetic interaction. 

On the border of ferromagnetism, magnetic pairing in the spin-triplet state has the dis- 
advantage that only the exchange of magnetic fluctuations polarised along the direction of 
the interacting spins, i.e., longitudinal fluctuations, contribute to the quasiparticle interac- 
tions. For a spin-rotationally invariant system, both longitudinal and transverse fluctuations 
contribute to pairing only for a spin-singlet state. 

Another disadvantage of being on the border of ferromagnetism is that for otherwise sim- 
ilar conditions the suppression of Tc due to the self interaction arising from the exchange of 
magnetic fluctuations is stronger than in the corresponding case on the border of antiferro- 
magnetism. This disadvantage can be mitigated in systems with strong magnetic anisotropy 
in that the effect of the transverse magnetic fluctuations on the self interaction would be 
suppressed while the strength of the pairing interaction arising from the longitudinal mag- 
netic fluctuations need not be reduced. This may apply in systems with strong spin-orbit 
interactions or in the spin-polarised state close to the border of ferromagnetism. 

These arguments have stimulated a new search for evidence of superconductivity on 
the border of itinerant electron ferromagnetism in cases where spin anisotropy is expected 
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to be pronounced, such as UGe2- This search has proved fruitful because it led to the first 
observation of the coexistence of superconductivity and itinerant electron ferromagnetism 
in UGe2 |§] and shortly thereafter in ZrZn2 and URhGe [^. 

The prediction of the simple model presented in ref. [0] that magnetic pairing is more 
robust in the quasi-2D square lattice than in the cubic structure seems to have been borne out 
by recent experiments. Namely, one finds an order of magnitude increase in the maximum 
Tc and in the range in pressure where superconductivity is observed on the border of metallic 
antiferromagnetism when the simple cubic lattice of Celn^ [§| is stretched along one principal 
axis by the insertion of non- magnetic layers to form the tetragonal compounds CeMIn^ 
(M is Co, Rh or Jr). 

These systems, albeit quite anisotropic, would not normally be considered to be quasi- 
two dimensional and it is not clear at first sight that the model calculations carried out in 
ref. are directly relevant. The purpose of this paper is to show that for comparable model 
parameters, the robustness of magnetic pairing increases gradually as one goes from a cubic 
structure to a more and more anisotropic structure on the border of metallic antiferromag- 
netism and ferromagnetism. This behaviour of the mean field transition temperature is in 
stark contrast to that of the "one- loop" fluctuations corrections to T^. The latter corrections 
typically depend logarithmically on the degree of anisotropy and would be expected to be 
negligible for materials such as CeMIn^. 

We do not expect some of the results of ref. to be generic properties of the magnetic 
interaction model. We have already stressed that even in simple cases, the robustness of 
magnetic pairing can be very sensitive to certain details of the lattice and electronic structure. 
Even in the single band problem, many such structures have not yet been extensively studied 
theoretically. Furthermore, we expect that the range of possibilities to be greatly expanded 
in the presence of more than one partially filled electronic band. 

Most known materials on the border of magnetism crystallise in other than simple cubic 
or tetragonal structure and have more than one band crossing the Fermi level. For these 
more complex systems, one would not expect the model of ref. to be directly relevant. 
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For example, the observation of spin-triplet rather than spin-singlet d-wave pairing in some 
multi-band materials with strongly enhanced antiferromagnetic spin fluctuations, such as 
UPts and Sr2Ru04, may not be inconsistent with the idea of magnetic pairing. A detailed 
study of magnetic pairing in multi-band systems for a range of crystal structures would shed 
light on the possible forms of superconductivity and the conditions most favorable for their 
observation. 

The simple model calculations suggest that anisotropic forms of superconductivity should 
be a generic property of systems on the border of metallic magnetism. It may seem surprising 
therefore that there are still so few observations of this phenomenon. In many cases, the 
multiplicity of bands and, for example, magnetic fluctuations in the non-bipartite lattice 
may weaken magnetic pairing to such an extent that quenched disorder may completely 
suppress superconductivity. An illustration of this point is the dramatic collapse of the spin- 
triplet superconducting transition temperature in Sr2RuOi in the presence of Al impurity 
concentrations as low as 0.1%. 

At flrst sight, the magnetic interaction model is mathematically analogous to the con- 
ventional electron-phonon problem with the generalised magnetic susceptibility playing the 
role of the phonon propagator. One would therefore expect that a simple analytic expression 
similar to that proposed by McMillan could be used to represent approximately the cal- 
culated numerically via the Eliashberg equations. Our attempts in this direction have not, 
however, proved successful A recent study suggests that there may be a fundamental 



reason for inapplicability of the McMillan-style expression for |TT| . On the border of long- 
range magnetic order, the incoherent part of the electron Green function, which is ignored in 
the simplest treatment, plays a major role in the formation of the Cooper pair condensate. 
The traditional picture in which superconductivity arose from pairing of well deflned (weakly 
damped) quasiparticles appears inadequate on the border of metallic magnetism even in the 
mean fleld Eliashberg treatment. 

We note that in our model the coupling of the quasiparticles to the magnetic fluctuations 
is a phenomenological constant to be inferred from normal state properties that formally 



includes that part of the vertex correction which is local in space and time. Calculations 
have shown that the neglect of vertex corrections that are non-local in space and time is 



justified at least in some cases of physical interest . When the magnetic correlation length 
becomes sufficiently large, however, these neglected non-local vertex corrections (including 
superconducting phase fluctuations) may become important. Their effect on and the 
normal state properties are as yet incompletely understood. 

II. MODEL 

We consider quasiparticles in a simple tetragonal lattice described by a dispersion relation 

ep = -2t(cos(p^) + cos{py) + atCos{pz)) 

- 4t'(cos(pa;) cos{py) + at cos(p^) cos(p^) + at cos(py) cos{pz)) (2.1) 

with hopping matrix elements t and t'. at represents the electronic structure anisotropy 
along the z direction, at = corresponds to the quasi-2D limit while at = 1 corresponds to 
the 3D cubic lattice. For simplicity, we measure all lengths in units of the respective lattice 
spacing. In order to reduce the number of independent parameters, we take t' = 0A5t and 
a band flUing factor n = 1.1 as in our earlier work. 

The effective interaction between quasiparticles is assumed to be isotropic in spin space 
and is deflned in terms of the coupling constant g and the generalised magnetic susceptibility 
which is assumed to have a simple analytical form consistent with the symmetry of the lattice. 

(2-2) 

where k and Kq are the correlation wavevectors or inverse correlation lengths in units of the 
lattice spacing in the basal plane, with and without strong magnetic correlations, respec- 
tively. Let 

g| = (4 + 2am) ± 2(cos(ga;) + cos{qy) + am cos{q^)) (2.3) 
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where am parameterises the magnetic anisotropy. = corresponds to quasi-2D magnetic 
correlations and a^, = 1 corresponds to 3D magnetic correlations. 

In the case of a nearly ferromagnetic metal the parameters and //(g) in Eq. ( |2.2| ) are 
defined as 

f = t (2.4) 
r^iq) = Tsfq- (2.5) 

where Tsj is a characteristic spin fluctuation temperature. Note that our definition of T,/ 
may differ from the characteristic spin fluctuation temperature scales used by other authors. 

In the case of a nearly antiferromagnetic metal, the parameters (f and 77(g) in Eq. ( p.2|) 
are deflned as 

f = ql (2.6) 

r/(g) = T,jq_ (2.7) 

As in our previous work |]T],0], the band structure and generalized magnetic susceptibility 
are modeled independently. This choice may be inconsistent when all of the contributions to 
x(q, uo) come from the chosen band. However, it allows us, in principle, to deal with the case 
where there are other important contributions to the generalized magnetic susceptibility. It 
has been argued that the latter case is of relevance to the ruthenates [|^, and most likely 
the heavy-fermion systems. 

A complete description of the model, the Eliashberg equations for the superconducting 
transition temperature and their method of solution can be found in the appendix. 

We note that the model is fully deflned by the phenomenological parameters describing 
the electronic structure ep, the generalised magnetic susceptibility ^(q, ^) and the interaction 
vertex g. In principle, these parameters can be estimated from experimental studies of the 
normal state. In particular, the resistivity can be used to estimate the dimensionless coupling 
parameter g^x^l'^ the value of which is between 10 and 20 for the simplest RPA model for 
the magnetic interaction. 



III. RESULTS 



A. Solution of the Eliashberg Equations for Tc 

The dimensionless parameters at our disposal are g^x^jt^ Tgf/t, kq and k. For comparison 
with the resuhs of our earher work [0,0, we take Tgf = |t and Kq = 12. In 2D, this 
Tsf corresponds to about 1000°/^ for a bandwidth of 1 eV while our choice of of is a 
representative value. 

The results of our numerical calculations of the mean field critical temperature as a 
function of the electronic and magnetic anisotropy parameters at and am, respectively, are 
shown in Figs 1 and 2 for representative values of the parameters and g'^Xo/t- Figures 
la-c illustrate the results for a nearly anti-ferromagnetic metal and figures 2a-c for a nearly 
ferromagnetic metal. Note that our previous calculations correspond to the quasi-2D case 
= ctm = and to 3D case at = a^a = 1- 

A glance at figure 1 reveals a clear pattern in the variation of Tc with both anisotropy 
parameters at and a^. We notice that T^. increases gradually and monotonically as the 
system becomes more and more anisotropic in either the electronic structure or in the mag- 
netic interaction. In going from 3D to quasi-2D, Tc/T^j is found to increase by up to an 
order of magnitude for otherwise fixed parameters of the model. The increase becomes least 
pronounced as for small and large (7^X0/^- 

The behavior in the nearly ferromagnetic case, figure 2, though broadly similar to that 
of the nearly antiferromagnetic metal, shows some interesting differences. In some cases, the 
minimum Tc occurs for 3D electronic structure, but quasi-2D magnetic interaction. Also, 
in all cases considered the maximum Tc is obtained for a quasi-2D electronic structure and 
strongly anisotropic, but not 2D magnetic interactions. 
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B. Mass renormalisation and interaction parameter 



In order to make a comparison with the corresponding electron-phonon problem it is 
instructive to define a mass renormahzation parameter \z and interaction parameter Aa- 
We define 

< -1- >FS{p) 

^ I-^ ^ < i^'i'V^jP P'- ^')//(p)//(pO >i-.s(p.pO .32) 

^ < ?72(p) >^5(p) 

where 

Vz{q,u;)=g\{q,u;) (3.3) 

and 

Vp(q,u;)^-^x(ci,oo) (3.4) 
7;(p) = sin(p^) (3.5) 

for p-wave spin triplet pairing (A = p) while 

Vd{ci,uj) = g\{q,uj) (3.6) 
77(p) = cos(p^) - cos{py) (3.7) 

in the case of d-wave spin-singlet pairing (A = d). The Fermi surface averages are given by 

• • • 5(ep - /x) (3.8) 

/d'^p d'^p' 
(2^(2^'"'^^^P~^^^^'P'~^^ ^^'^^ 

In practice, we compute the Fermi surface average with a discrete set of momenta on a cubic 
or tretragonal lattice and we replace the delta function by a finite temperature expression 

■Slfp - — ^/p(l - /p) (3.11) 



where /p is the Fermi function. Note that ^fp{l — /p) — > S{ep — ji) as T ^ 0. We have used 
T = O.lt and = 128'^ in all of our calculations. The finite temperature effectively means 
that van Hove singularities will be smeared out. 

Note that the Fermi surface average that appears in A^, Eq. plays a role similar to 
that of a^F{ijj)/uj in the case of phonon mediated superconductivity. From the definitions 
of the parameters \za ( F^ , (|3^ ) and our model for f^) Eq- (|2.2|) , we see that 
Az,A are directly proportional to the dimensionless factor g'^Xoi'^o/t- Thus we will consider 
the quantities 

Az,A = XzA/{9\o^l/t) (3-12) 

which are functions only of n, t' /t and k^. 

In figures 3 and 4 we show A^, A^ and the ratio Aa/A^ for a representative value of k"^ in 
the case of a nearly antiferromagnetic metal and nearly ferromagnetic metal, respectively. 

The trends in both cases are the same. A| and A^ are seen to increase gradually and 
monotonically in going from 3D to quasi-2D. However, A^ grows faster than A^ so the 
ratio Aa/A^ also increases in going from 3D to quasi-2D. This qualitative trend in the 
ratio is consistent with the behavior of Tc obtained from the numerical solution of the 
Eliashberg equations. In the ferromagnetic case, however, it fails to reproduce the fact that 
the minimum Tc is not necessary for a fully 3D system and that the maximum Tc is obtained 
for strongly anisotropic yet not quasi-2D systems. 

IV. DISCUSSION 

The results of the calculations for both the nearly ferromagnetic and nearly antiferro- 
magnetic metals show that the robustness of magnetic pairing increases gradually as one 
goes from a cubic to a more and more anisotropic structure with parameters other than 
and at left unchanged. These results are consistent with our previous findings [0 and with 
the calculations for am = at = and am = ctt = 1 presented in ref. [|I^]. In an earlier 
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study, Nakamura et al. found that Tc could increase by up to a factor of three in 

going from 3D to 2D for their choice of model parameters. The effect of anisotropy on Tc 
for nearly ferromagnetic and nearly antiferromagnetic metals is qualitatively similar. This 
phenomenon arises from the increase with growing anisotropy of the density of states of both 
the quasiparticles and of the magnetic fluctuations that mediate the quasiparticle interac- 
tion. This effect could be further enhanced in the case of a nearly antiferromagnetic metal 
by the change in the pattern of the oscillations of the magnetic interaction. 

It can be seen from figure 5 that the strength of the interaction in the repulsive sites 
outside of the nodal plane of the dx2_y2 state gets reduced while crucially the attraction in 
the basal plane gets enhanced as one goes from the cubic to a more and more anisotropic 
tetragonal lattice. This enhancement is the consequence of the increase of the phase space 
of soft magnetic fluctuations as one goes from a cubic to a quasi-two dimensional structure. 
Since our model potential varies smoothly with the tetragonal distortion, parameterised by 
am in figures 1 and 2, it is clear that these effects occur gradually with increasing separation 
between the basal planes. 

The calculations assume that the maximum magnetic response for a nearly antiferro- 
magnetic metal occurs at the commensurate wavevector defined hj Qx = Qy = Tc/a and 
Qz = 7r/c, where a and c are the lattice constants in the basal plane and along the tetrag- 
onal axis respectively, reintroduced here for clarity. The oscillations in the magnetic inter- 
action potential along the tetragonal axis obviously depend on the value of Qz- However, 
the enhancement of the attraction in the basal plane and the reduction of the interaction 
elsewhere as one goes from a cubic to a more and more anisotropic lattice do not depend on 
the particular value of Qz- Therefore, we expect the qualitative conclusions of this paper to 
be independent of Qz- 

The robustness of the pairing is further enhanced by the gradual change in the electronic 
band from a 3D to a quasi-2D form (see Eq. The reduced hopping along the distortion 

axis, parameterised by at in figures 1 and 2, implies a reduced electronic bandwidth and 
hence increased density of electronic states. Our calculations show that this too leads to a 
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gradual increase in Tc with increasing distortion of the lattice. 

In a nearly ferromagnetic metal, one again benefits from the reduction of the electronic 
band width and the increase of the interaction in the basal plane as one goes from a cubic 
to a tetragonal lattice (see figure 6). However, the suppression of the interaction between 
the basal planes has a less dramatic effect on the border of ferromagnetism than antiferro- 
magnetism because in the latter case one suppresses key repulsive regions of the interaction 
(figure 5). 

These simple arguments explain how the pairing effects of the interaction are strength- 
ened by a tetragonal distortion in our model. However, the same effects also contribute 
to an enhanced self-interaction which acts to suppress T^. The relative importance of the 
pair forming and pair breaking effects of the magnetic interaction cannot be inferred by 
the above physical picture alone. The numerical calculations show that for most cases con- 
sidered here the pair forming effects dominate. The balance is particularly delicate on the 
border of ferromagnetism where the suppression of T^. brought about by the self-interaction 
is pronounced. A physical interpretation of this suppression of T^. is given in ref. |]TT|. The 
same interpretation may explain, for example, why the maximum of Tc/T^f in the nearly 
ferromagnetic case is for a strongly anisotropic yet not quasi-2D pairing potential (figure 2). 

A most striking manifestation of the interplay between the pair-forming and pair- 
breaking tendency of the magnetic interaction is the breakdown of the McMillan-style ex- 
pression for Tc in terms of the parameters A a and \z (see Eqs. (p.l|J3.2| ). This was noted 
in ref. and has been interpreted in ref. in terms of the important role played by 
the incoherent part of the Green function which is ignored in the simplest treatments, but 
is included in the present and earlier work where the full momentum and frequency 
dependence of the self energy is taken into account. 
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V. OUTLOOK 



The calculations show that the lattice anisotropy may increase the robustness of magnetic 
pairing in the mean-field approximation. Superconducting phase fiuctuations which are not 
included in this approximation may be expected to suppress Tc in the 2D limit. Therefore, 
in practice, one would think that the most favorable case for magnetic pairing is that of 
strong but not extreme anisotropy. 

As noted in the introduction and in the previous two sections, the robustness of magnetic 
pairing can be very sensitive to certain details of the magnetic interaction and electronic 
structure. Therefore, one should exercise caution in making quantitative comparisons be- 
tween the results of our calculations and experiment. For instance, one would expect all of 
the parameters of the model (not solely and at) to change simultaneously with increas- 
ing lattice anisotropy. The changes brought about in going from a cubic to a tetragonal 
lattice may even be much more complex than considered here. In particular, the number of 
partially filled bands may itself change. As also mentioned in the introduction, this could 
have in some cases even more dramatic consequences on superconductivity than the effects 
taken into account in our simple one-band model. 

The theoretical framework developed for systems on the border of magnetism can be 
translated to describe systems on the border of other types of instabilities, such as charge 
density wave or ferroelectric instabilities. The above given phase space argument to explain 
the increased robustness of magnetic pairing with increasing lattice anisotropy should carry 
over in part to these other pairing mechanisms, at least at the one-loop mean-field level (see. 



e.g., ref. |[T5[] ). 

While some understanding of the properties of the magnetic interaction model has been 
gained over the last few years (e.g., the conditions for robust pairing of electrons), there are 
many cases where the predictions of the model have not been worked out. Of particular 
importance is the role of the multiplicity of partially filled bands which may be expected to 
be the key to understanding exotic superconductivity observed in nearly magnetic materials 
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such as f/Pts and Sr2RuOi. 
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VII. APPENDIX 

We consider quasiparticles on a cubic or tetragonal lattice. We assume that the dominant 
scattering mechanism is of magnetic origin and postulate the following low-energy effective 
action for the quasiparticles: 

Seff = Y. dripl^ (r) [dr + Cp - /x) ^p,^ (r ) 

-^E //rfr^'rfr'x(q,r-r')s(q,r)-s(-q,r') (7.1) 

where N is the number of allowed wavevectors in the Brillouin Zone and the spin density 
s(q, r) is given by 

E ^P+q,"(^)^":7'^P,7(^) (''•2) 
p,o,7 

where a denotes the three Pauli matrices. The quasiparticle dispersion relation ep is defined 
in Eq. ( p.l|) , /i denotes the chemical potential, /? the inverse temperature, g'^ the coupling 
constant and ipp u V'p.o- are Grassmann variables. In the following we shall measure 
temperatures, frequencies and energies in the same units. 

The retarded generalized magnetic susceptibility x{s\.i ^^^^ defines the effective inter- 
action, Eq. (|7?1|), is defined in Eq. (p.2[). 

The spin-fluctuation propagator on the imaginary axis, x(q, iv^i) is related to the imag- 
inary part of the response function Jmx(q, cj), Eq. (|2.2|) , via the spectral representation 
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Xq,«i^„ = -/ ■■ (7-3 

To get x(q5 i^n) to decay as as z/„ oo, as it should, we introduce a cutoff c^o and take 
/mx(q, tu) = for > wq. A natural choice for the cutoff is = ri{q)Kl. We have checked 
that our results for the critical temperature are not sensitive to the particular choice of ujq 
used. 

The Eliashberg equations for the critical temperature Tc in the Matsubara representation 
reduce, for the effective action Eq. ([7.1[), to 

T 

S(p, iujn) = g^j^ X(P - k, iujn - i^n)Gik, iVln) (7.4) 



G{p,iuJn) = -. J ^ —, : r (7.5) 

lUJn - (Cp - A*) - S(p, lUJn) 



9l 

3 

i2 



T 

]^ X(P - k, ioJn - ifi„) |G(k, ifi„) p<l>(k, zfi. 



A(T) = 1 — >T = T^ (7.6) 

where E(p,ia;„) is the quasiparticle self-energy, G{p,iujn) the one-particle Green's function 
and $(p, -itUn) the anomalous self-energy. The chemical potential is adjusted to give an 
electron density of n = 1.1, and is the total number of allowed wavevectors in the 
Brillouin Zone. In Eq. (|7.6|) , the prefactor g'^/3 is for triplet pairing while the prefactor —g"^ 
is appropriate for singlet pairing. Only the longitudinal spin-fluctuation mode contributes 
to the pairing amplitude in the triplet channel. Both transverse and longitudinal spin- 
fluctuation modes contribute to the pairing amplitude in the singlet channel. All three 
modes contribute to the quasiparticle self-energy. 

The momentum convolutions in Eqs. ( [7.4| , |7.6| ) are carried out with a Fast Fourier Trans- 
form algorithm on a 48 x 48 x 48 lattice. The frequency sums in both the self-energy 
and linearized gap equations are treated with the renormalization group technique of Pao 



and Bickers []TB[. We have kept between 8 and 16 Matsubara frequencies at each stage of 



the renormalization procedure, starting with an initial temperature Tq = 0.6t and cutoff 
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~ sot. The renormalization group acceleration technique restricts one to a discrete set of 
temperatures Tq > Ti > T2 . . .. The critical temperature at which A(T) = 1 in Eq. ( [r.6| ) is 
determined by linear interpolation. 



16 



REFERENCES 

[1] P. Monthoux and G.G. Lonzarich, Phys. Rev. B 59, 14598 (1999). 

[2] P. Monthoux and G.G. Lonzarich, Phys. Rev. B 63, 054529 (2001). 

[3] T.Moriya, Y. Takahashi and K. Ueda, J. Phys. Soc. Jpn. 52, 2905 (1990). P. Monthoux, 
A.V. Balatsky and D. Pines, Phys. Rev. Lett. 67, 3448 (1991). 

[4] E.R. Dobbs, "Helium Three", Oxford University Press (2001). 

[5] S.S. Saxena et al.. Nature 406, 587 (2000). 

[6] C. Pfleiderer et al. Nature 412, 58 (2001). 

[7] A. Aoki et al. Nature 413, 613 (2001). 

[8] I. R. Walker, P.M. Grosche, D. M. Freye, and G.G. Lonzarich, Physica C 282, 303 
(1997). 

[9] C. Petrovic et al., J. Phys. C 13, L337 (2001). 

[10] T. Kuwabara and M. Ogata, Phys. Rev. Lett. 85, 4586 (2000). 

[11] A. Abanov, A.V. Chubukov, and A.M. Finkel'stein, Europhys. Lett. 54, 488 (2001). 

[12] P. Monthoux, Phys. Rev. B 55, 15261 (1997). A.V. Chubukov, P. Monthoux and D.R. 
Morr, Phys. Rev. B 56, 7789 (1997). P. Monthoux, Phil. Mag. B 79, 15 (1999). S. 
Nakamura, J.Phys. Soc. Jpn. 69, 178 (2000). N.E. Bickers and S.R. White, Phys. Rev. 
B 43, 8044 (1991). 

[13] R. Arita, K. Kuroki, and H. Aoki, Phys. Rev. B 60, 14585 (1999) ; J. Phys. Soc. Jpn. 
69, 1181 (2000). 

[14] S. Nakamura, T. Moriya and K. Ueda, J. Phys. Soc. Jpn. 65, 4026 (1996). 
[15] I.G. Khalil, M. Teter and N.W. Ashcroft, Phys. Rev. B 65, 195309 (2002). 



17 



[16] C.-H. Pao and N.E. Bickers, Phys. Rev. B 49, 1586 (1994). 



18 



FIGURES 

Nearly Antiferromagnetic: = 0.25 ; Xo /t = 5 
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FIG. 1. Eliashberg Tc/Tgf for nearly antiferromagnetic systems as a function of the electronic 
anisotropy parameter at and the magnetic anisotropy parameter for representative values of 
the correlation wavevector and coupling constant g^Xo/'t- (a) = 0.25, g^Xa/'t = 5. (b) 

= 0.50, g^Xo/t = 10. (c) = 1.00, g'^xo/t = 10. = = corresponds to the 2D limit 
while at = am = 1 corresponds to an isotropic 3D system. 
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Nearly Ferromagnetic: K = 0.25 ; g Xo = 5 



(a) 




Nearly Ferromagnetic: = 0.50 ; g Xo = 10 
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Nearly Ferromagnetic: = 1.00 ; g Xo = 10 
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FIG. 2. Eliashberg Tc/Tgf for nearly ferromagnetic systems as a function of the electronic 
anisotropy parameter at and the magnetic anisotropy parameter am for representative values of 
the correlation wavevector and coupling constant ff^Xo/*- (a) = 0.25, g^Xo/t = 5. (b) 



0.50, ff^^o/i = 10. (c) = 1.00, g\o/t = 10. 



Or 



corresponds to the 2D limit 



while at = am = 1 corresponds to an isotropic 3D system. 
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Nearly Antiferromagnetic: K = 0.25 
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Nearly Antiferromagnetic: K. = 0.25 
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FIG. 3. Interaction parameters (a) A^, (b) and ratio (c) \dl^z for nearly antiferromagnetic 
metals for a representative value of = 0.25 



Nearly Ferromagnetic: K = 0.25 
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Nearly Ferromagnetic: K = 0.25 
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Nearly Ferromagnetic: K = 0.25 

(c) 
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FIG. 5. The magnetic potential seen by a quasiparticlc in a spin-singlet dx2-y2 Cooper pair 
state given that the other quasiparticle is at the origin (marked by a cross) . The figure depicts the 
evolution of the potential as one goes from a cubic to a tetragonal lattice by varying the parameter 
am- Closed circles denote repulsive sites and open circles attractive ones. The size of the circle is 
a measure of the strength of the interaction. The nodal plane of the d^2_y2 state are represented 
by the shaded region. 
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FIG. 6. The magnetic potential seen by a quasiparticlc in a spin-triplet px Cooper pair state 
given that the other quasiparticle is at the origin (marked by a cross). The figure depicts the 
evolution of the potential as one goes from a cubic to a tetragonal lattice by varying the parameter 
am- Open circles denote attractive sites. The size of the circle is a measure of the strength of the 
interaction. The nodal plane of the Px state is represented by the shaded region. 
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